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Introduction 

Fast  fracture  induces  rapid  particle  motions  in  the  cracking  body. 

From  the  point  of  view  of  fracture  mechanics  such  dynamic  effects  are  of 
particular  interest  in  the  immediate  vicinity  of  the  propagating  crack  tip, 
where  the  fields  of  stress  and  deformation  are  critical  for  either  con¬ 
tinued  fracture  or  for  crack  arrest. 

A  substantial  body  of  literature  has  dealt  with  dynamic  effects  on 
essentially  brittle  fracture,  within  the  context  of  linear  elastic  fracture 
mechanics.  We  mention  review  articles  by  Achenbach  (1972),  Freund  (1975) 
and  Kanninen  (1978).  Dynamic  effects  in  the  presence  of  elastic-plastic 
constitutive  behavior  have,  however,  been  considered  in  only  a  few  studies. 

An  investigation  of  the  dynamic  near-tip  fields  in  an  elastic  perfectly- 
plastic  material  was  presented  by  Slepyan  (1976),  who  considered  both  the 
cases  of  anti-plane  and  in-plane  strains.  Dynamic  near-tip  effects  in  the 
presence  of  strain  hardening  were  investigated  by  Achenbach  and  Kanninen 
(1978)  on  the  basis  of  J ^  flow  theory.  These  authors  found  results  which 
are  very  similar  to  the  ones  obtained  by  Amazigo  and  Hutchinson  (1977)  for 
the  corresponding  quasi-static  problem.  As  shown  by  Achenbach,  Burgers  and 
Dunayevsky  (1979),  for  strain  hardening  the  governing  equations  are  elliptic 
when  the  crack-tip  speed  is  less  than  a  certain  critical  value.  The  usual 
separation-of-variables  asymptotic  analysis  can  then  be  carried  out,  which 
yields  singularities  of  the  general  type  r**(-l  <  p  <  0)  for  the  stresses 
and  the  strains.  As  the  crack-tip  speed  increases  (or  alternatively  as  the 
strain-hardening  curve  becomes  flatter)  the  nature  of  the  governing  equa¬ 
tions  becomes,  however,  hyperbolic,  and  the  near-tip  fields  appear  to  change 


character.  Indeed  in  the  limit  of  elastic  perfectly-plastic  behavior  the 
stresses  become  bounded  and  only  some  strains  display  singularities,  as 
shown  in  the  present  paper,  and  earlier  by  Slepyan  (1976). 

In  this  paper  dynamic  effects  on  near- tip  fields  are  investigated  for 
elastic  perfectly-plastic  constitutive  behavior.  The  approach  is  different 
from  the  one  employed  by  Slepyan  (1976),  but  the  results  are  identical  for 
the  Mode-Ill  case.  For  Mode-I  the  results  show  some  differences.  In  the 
present  paper  the  stresses  and  the  strains  have  been  expressed  in  explicit 
form,  and  they  include  higher  order  terms. 

As  the  crack-tip  speed  decreases  the  expressions  for  the  stresses  reduce 
to  the  ones  for  the  corresponding  quasi-static  problem,  as  might  be  expected  on 
the  basis  of  intuitive  arguments.  This  is,  however,  not  true  for  the  strains, 
which  become  unbounded  in  the  limit  of  vanishing  crack-tip  speed.  A  signi¬ 
ficant  difference  between  the  dynamic  and  quasi-static  solutions  is  that  the 
dynamic  solution  does  not  describe  elastic  unloading  behind  the  crack  tip. 

These  anomalies  suggest  that  the  transition  from  dynamic  to  quasi-static 
conditions  is  non-uniform.  It  is  speculated  that  the  transition  from  the 
dynamic  to  the  quasi-static  solution  with  decreasing  crack- tip  speed  is 
effected  because  the  dynamic  solution  is  asymptotically  valid  in  a  small 
zone,  which  shrinks  on  the  crack  tip  in  the  limit  of  vanishing  crack-tip 
speed . 

1.  Governing  Equations 

Both  a  stationary  coordinate  system  with  axes  denoted  by  x^ ,  and  a 
moving  coordinate  system  with  axes  denoted  by  (x,y,z)  are  considered.  The 
moving  coordinate  system  has  its  origin  at  the  propagating  crack  tip.  The 
geometry  is  shown  in  Fig.  1.  In  this  section  the  equations  governing  the 


motions  of  an  elastic  perfectly-plastic  material  are  stated  in  the  stationary 
coordinate  system.  In  the  next  sections  these  equations  are  simplified  for 
anti-plane  strain  and  plane  strain,  for  the  special  case  of  "steady-state” 
fields  of  stress  and  deformation  relative  to  the  moving  crack  tip. 

Relative  to  the  stationary  coordinate  system  the  equations  of  motion 


are 


d.tf,  . 

J 


,2 

P  o^u 


t  i 


(1.1) 


In  the  zone  of  plastic  deformation  the  stresses  are  assumed  to  satisfy 
the  Tresca  yield  condition,  which  states  that 

|tI  =  k  (1.2) 

where  lTlmax  Is  the  maximum  shear  stress,  and  k  is  the  yield  stress  in  pure 
shear.  For  an  elastic  perfectly  plastic  solid  the  total  strain  rates  are 
defined  by 

5  e. .  =  S  e?1  +  a  s?1  (1.3) 

t  ij  t  ij  t  ij  ^  ' 


Here  the  elastic  strain-rates  are  defined  by 

Vij  =  2^  (Vij  ‘  ~  Vkk6ij) 

while  the  plastic  strain  rates  are 

.P1 


t  ij 


\  s 


ij 


(1.4) 


(1.5) 


Equations  (1.3)  -  (1.5)  are  the  Prandtl-Reuss  flow  equations 


In  (1.4)  y,  and  v  are  Che  elastic  shear  modulus  and  Poisson's  ratio. 


respectively,  while  in  (1.5)  s„  defines  the  stress  deviator,  and  \  is  a 
non-negative  proportionality  factor,  which  may  vary  in  space  and  time. 
Equations  (1.1)  -  (1.5)  should  be  supplemented  by  appropriate  boundary 
conditions . 

The  speed  of  the  crack  tip  is  v  =  di(t)/dt,  where  2(t)  is  the  rnono- 
tically  increasing  crack  length,  and  where  v(t)  and  dv/dt  are  continuous. 
The  moving  and  the  stationary  coordinate  system  are  related  by 

x  =  xL  -  2 (t)  (1.6) 

By  means  of  (1.6)  the  material  time  derivatives  are  transformed  to  the 
moving  coordinate  system  by  the  relations 


V  >  =  IF  - 


(1.7) 


-t(  )  =  7"T  '  *(t)^  '  2v(t)^  +  v2(t)7T 

Ot  -X 


(i.8) 


For  the  case  that  v(t)  approaches  a  constant  value  v  as  t  increases, 
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it  may  be  assumed  that  steady-state  fields  of  stress  and  deformation 
are  established  relative  to  the  coordinate  system  moving  with  the  crack 
tip.  This  assumption  implies  that  the  time  derivatives  on  the  right- 
hand  sides  of  (1.7)  and  (1.8)  may  be  neglected,  which  leads  to  the  usual 
results  that 


a  2  2  V 

>  C  \  ~  _  v  \  ~  _s_ 

V  '  »  *x  ’  V  ’  -  2 

^x 


(1.9a,b) 


In  the  sequel  the  speed  of  the  crack  tip  will  appear  in  the  "Mach 


number"  M  ,  which  is  defined  as 


(1.10) 


M  =  v  J{v.l  o) 

where  (jjk/o)  defines  the  velocity  of  shear  waves  in  an  elastic  solid 
with  shear  modulus  y, . 

2 .  Crack- Propagation  in  Anti-Plane  Strain 
In  the  moving  coordinate  system  (x,v,z),  steady-state  motion  in  anti¬ 
plane  strain  is  defined  by  a  displacement  w(x,y)  in  the  z-direction.  3v 
using  (1.9b),  the  equation  of  motion  (1.1)  then  reduces  to 


2 

+3  o'  -ov  5  w  =  0 
x  xz  y  yz  =  xx 


The  Tresca  yield  condition  (1.2)  may  be  expressed  in  the  form 


and 


the  Prandtl-Reuss  flow  law  (1.3)  yields  by  using  (1.9a) 

5  w  =  -  3  c  -2  \  g 

XX  ;*  X  xz  xz 


(2.1) 


(2.2) 


(2.3) 


xy  n  x  yz  yz 

It  is  sufficient  to  consider  the  solution  in  the  half-plane  y  s  0  . 
The  boundary  conditions  at  y  =  0  are 


X  <  0  : 

a  =0 
yz 

x  2  0  : 

c  >0 
yz 

,  w  =  0 

The  yield  condition  (2.2)  is  identically  satisfied  by 


a  =  -k  sinw  ,  c  =  k  cosui 
xz  yz 


By  introducing  the  strain  component 


w  =  p  w 
x  x 


(2.4) 


(2.5) 

(2 . 6a , b ) 


(2. 7a, b) 


(2.8) 


together  with  (2.7a,b)  in  (2.1),  and  eliminating  2\  from  (2.3)  and 
(2.4),  we  obtain 


2 

cosiu  3  jj  +  siniu  :  *  +  M  j1  :  w  =  0  (2.9) 

x  y  k  x  x 

and 

cos«  )  w  +  sinw  )w  r-;ji  =  0  (2. 10) 

xx  y  x  u  x 

where  M  is  defined  by  (1.10). 

Equations  (2.9)  and  (2.10)  constitute  a  hyperbolic  system  of  equations. 
A  brief  study  of  (2.9)  and  (2.10),  including  a  discussion  of  the  charac¬ 
teristic  curves  and  the  corresponding  Riemann  invariants  has  been  pre¬ 
sented  by  Achenbach,  Burgers  and  Dunayevsky  (1979).  These  authors  showed 
that  a  solution  of  (2.9)-(2.1Q)  satisfying  the  boundary  conditions  (2.5) 
and  (2.6a,b)  can  be  constructed  by  using  only  simple  and  uniform  fields. 

In  the  present  paper  we  start  off  with  so-called  "Simple  wave"  solutions, 
which  are  solutions  for  which  the  dependent  variables  are  functions  of 
each  other,  a  brief  discussion  of  simple  wave  solutions  is  presented  in 


Appendix  A. 

It  is  easily  verified  that  (2.9)  and  (2.10)  are  of  the  general  form 
given  by  (A.l),  where 


(2 . 11a  ,b ) 


(2.12) 


The  characteristic  directions  follow  from  (A. 7)  as 


for  v  :  w  =  -  (k/iiM)  *  +  w  (2.18) 

x  x 

where  w  and  w  are  constants  of  integration, 
xx 

Equation  (2.17)  leads  to  a  negative  value  of  V.  Hence  we  use  only 
the  relation  given  by  (2.18).  Substitution  of  (2.18)  into  either  (2.9) 
or  (2.10)  gives 

(cos*  -  M)  3  *  +  sin*  3  *  *  0  (2.19) 

x  y 

Along  the  family  of  characteristic  curves  defined  by  V  ,  see  (2.13), 

(2.19)  reduces  to  d*  =  0  ,  i.e.,  *  =  constant.  Consequently  the  re¬ 
lations  given  by  (2.13)  for  the  V  family  can  be  integrated  to  yield 


y  -  - —  x 

J  cos*  -  M 


where  9  is  the  polar  angle  defined  by 
tan9  =  y/x  , 


(2.20) 


(2.21) 


and  t)(9)  is  an  arbitrary  function.  However,  by  virtue  of  the  boundary 
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condition  (2.6b)  we  must  have  3(5)  a  0  .  Thus 


smu. 

y  - - rr  x 

J  cosuj  -  M 


Equations  (2.22)  and  (2.21)  imply  a  relation  between  uj  and  5 
which  can  be  solved  to  yield 


2  2  2  k 

cosuj  =  M  sin  5  +  (1-M  sin  3)  coso 


(2.23) 


The  stresses  then  follow  from  (2.7a,b)  as 


-  2  2^ 

=  -  k  .(1-M  sin  3)  -  M  coso.  sine 


(2.24) 


=  k  _ ( 1-M^sin^3)^  coso  +  M  sin^S] 


(2.25) 


Substitution  of  (2.23)  into  (2.18)  yields 


-lr  2  22  k 

wx  ~  ~  (k/^M)  cos  LM  sin  S  +  (1-M  sin  3;  coso j 


(2.26) 


where  we  have  taken  w^  =  0  to  satisfy  the  boundary  condition  (2.6b). 
From  (2.10)  we  obtain 


q  w  =  3  w  = 


k  cosuj 


o  uj  - : —  o  « 


y  x  "x  y  uM  sinuj  ~x  sinuj  x 


(2.27) 


where  (2.18)  has  also  been  used,  and  is  defined  as 
w^  =  Bw/3y 


(2.28) 


Equation  (2.27)  can  be  integrated  with  respect  to  x,  to  yield 


k  ^  1-M  r  2  2  2  ^, 

w  =  ~  {~tT7  2hl1-M  sin  3  -  (1-M  sin  0)1  cosS] 
y  n  '  iM 


1+M  r  2  2  2  k 

+  ~2yi  2nLl+M  sin  3  +  (1-M  sin  3)  cos6jj  +  j(y) 


(2.29) 


where  Ky)  is  an  as  yet  undetermined  function. 


MjT*L-**f — ft 


I 
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The  function  \(r,9)  can  be  obtained  from  (2  24)  -  (2.26)  and  either 
(2.3)  or  (2.4).  After  some  manipulation  we  find 


*  r 


2  2  k 

(1  -  M  sin  c;  > 


(2.30) 


which  satisfies  the  essential  condition  on  \. 

It  is  noted  that  the  expression  for  a  given  by  (2.25)  cannot  satis- 

y  z 

fy  the  condition  that  CyZ  =  0  at  8  =  n  ,  as  required  by  (2.5).  This 
suggests  that  (2.23)  -  (2.30)  are  valid  only  in  a  domain  0  s  j  s  8  , 

k 

while  another  solution  holds  for  6  £  S  £  tt  .  This  other  solution  cannot 
be  a  simple  wave  solution,  since  no  other  than  the  one  given  by  (2.23)  - 

k 

(2.30)  is  available.  A  suitable  solution  for  &  £  0  £  n  is,  however, 

given  by  a  uniform  field  of  the  form: 

jj  =  constant,  w^  =  constant  (2.31) 

Such  a  field  clearly  satisfies  (2.9)  and  (2.10).  An  appropriate  choice 
for  ui  is 

*  =  -t/2  ,  (2.32) 

which,  by  virtue  of  (2.7a,b),  leads  to  the  following  stresses  for 

* 

6  £  9  £  TT 


■  -  k  ■  v  0  (2.33a,b) 

Thus,  the  boundary  condition  (1.5)  is  satisfied.  It  immediately  follows 
from  (2.23)  and  (2.32)  that  the  angle  9  must  satisfy  the  equation 

M  sin20*  +  (1  -  M2sin29*)^  cosG*  =  0  (2.34) 

Hence 

9*  =  -  tan'1  (1/M)  (2.35) 

Thus  0  varies  from  0  =  90°  at  M  =  0  to  9  *  135°  at  M  =  1  . 


Substitution  of  (2.35)  in  (2.26)  and  (2.29),  and  the  requirement 
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that  the  strains  should  be  continuous  at  8  =  G  leads  to  the  following 
results  for  6  £  6  £  n  ; 

w  ■  -  rk/2  n  M  (2.36) 


(2.37) 


Finally,  an  approximate  expression  for  j(y)  can  be  found  from  the 


condition  that  w  should  be  bounded  for  y  =  0,  x  ~  r  ,  where  r  is  the 

y  P  P 

length  of  the  plastic  zone.  For  small  y  the  expression  for  w  given  by 
(2.29)  then  implies  that 


*(y)  =  ~  2n(r  /|y|)  -  ■—  2n(l-M)  -  k  2  +  1  (2.38) 

The  results  presented  in  this  Section  are  essentially  the  same  as 

those  found  by  Slepyan  (19/6)  who  used  an  approach  in  which  only  the  most 

singular  terms  ,  0(l/r),  were  retained  in  the  governing  equations. 

For  three  values  of  M  the  stress  components  a  ,  given  by  (2.24),  (2.33a), 

and  a  given  by  (2.25),  (2.33b),  have  been  plotted  in  Figs.  2  and  3,  res- 
yz 

pectively.  The  curves  show  a  substantial  dependence  of  the  stresses  on  M. 

For  M  =  0,  (2.24)  and  (2.25)  reduce  to  a  =  -k  sin6  and  a  =  k  cos&  , 

xz  yz 

respectively,  which  are  the  expressions  for  the  quasi-static  stresses 
derived  by  Chitaley  and  McClintock  (1971).  Thus,  as  one  would  perhaps 
expect  intuitively,  the  dynamic  stresses  reduce  to  the  quasi-static  ones  in 
the  limit  M  —  0.  This  is,  however, not  true  for  the  strains  as  discussed 
later. 

For  0  <  6  £  0  ,  the  strain  (k*/k)w^  ,  which  is  given  by  (2.26),  has 
been  plotted  in  Fig.  4.  Again,  a  very  noticeable  dependence  on  M  is 
observed . 
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It  follows  from  (2.29)  and  (2.38)  that  w^  is  singular  at  the  crack 

•k 

tip.  In  the  domain  0  i  5  s  6  we  find 


k  1-M  ,  .  ,  N 

w  ~ - —  In  (r/r  ) 

y  n  M  p 


(2.39) 


This  result  shows  that  w^  not  only  becomes  unbounded  as  r  -  0  ,  but  also 
as  M  -  0.  The  strain  component  w^  is  bounded  as  r  —  0,  but  it  becomes 


unbounded  in  the  limit  M  -  0.  Thus,  for  the  strains  there  is  no  uniform 


transition  fromdynamic  to  quasi-static  conditions. 

In  the  loading  zone  the  quasi-static  solution  is  of  the  form 

w  ~  In(r/r  )sinS  as  r  "*  0.  Behind  the  crack  tip  the  quasi-static  solution  * 
x  p 

shows  a  region  of  elastic  unloading.  As  a  material  point  passes  the  crack 

tip,  the  increments  of  w^  change  sign  ,  and  for  the  quasi-static  case  this 

can  only  be  achieved  by  elastic  unloading.  For  the  asymptotically  valid 

dynamic  solution  presented  in  this  paper  the  loading  zone  extends,  however, 

completely  around  the  crack  tip,  even  though  different  expressions  hold  for 
★  * 

0  £  6  £  0  and  9  s  6  s  n,  respectively.  It  is  speculated  that  the  length 

dimensions  of  the  zone  in  which  the  present  dynamic  results  are  valid  are  of 
2 

order  0(M  ).  Outside  of  this  zone  another  (dynamic)  solution  which  does 
describe  unloading  and  which  does  show  a  uniform  transition  to  the  quasi¬ 
static  solution  may  occur. As  M  —  0  the  dynamic  zone  then  shrinks  on  the 
crack  tip,  and  in  the  limit  M  =  0  the  zone  has  vanished  altogether,  and  the 
quasi-static  solution  describes  the  near-tip  strains. 

A  final  result  of  interest  is  the  crack-opening  angle.  This  angle, 
which  is  defined  by  tan^  a  -  \w^ |  at  9  »  tt  ,  follows  from  (2.36)  as 

tan  a  =  2  tan’1  (nk/2^M)  (2.40) 

As  M  —  0  we  find  a  ■  tt  ,  which  is  in  agreement  with  the  quasi-static  result. 


3 .  Crack- Propagation  in  Plane  Strain 


In  the  moving  coordinate  system  (x,y,z),  steady-state  motion  in  plane 

strain  is  defined  by  displacements  u(x,y)  and  v(x,y)  in  the  x-  and  y- 

cirections,  respectively.  We  assume  that  G  is  the  intermediate  stress, 

z 

i.e.,  a  £  g  <.  g  .  For  the  elastic  case  this  is  true,  for  the  elastic- 
x  z  y 

plastic  problem  it  must  be  checked  a  posteriori.  It  was  noted  by  Koiter 
(1953),  see  also  Gelring®  (1973,  p.  514),  that  the  Tresca  yield  condition 
then  implies  that 


^  -  gP^  =  0 

"t  z  z 


(3.1) 


el 


The  condition  of  plane  strain  consequently  reduces  to  e2  =  0  ,  which 
by  virtue  of  (1.4)  yields  the  result 

G  =  V(G  +  G  ) 

z  x  y 

We  now  introduce  the  following  new  variables 

•  I y  •  y 

<?  =  7  (a  +a  +a) 

3  x  y  z 

u  =  S  u 

X  X 

V  =  ?  V 
X  X 


(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.6) 


By  the  use  of  Eq.(3.2)  it  is  then  easily  checked  that  the  equations  of 
motion  in  the  moving  coordinate  system  may  be  expressed  in  the  form 


S  °  +  o / iT, T  £a  +  da  -  p  v2  *  u  =  0 

x  -  2(l+v)  x  y  xy  ®  x  x 


(3.7) 


-V-+  2(l+v)  V 


3  2 

d  a  +  ~  G  -  o  v  3  v 
X  xy  31  X  X 


(3.8) 


Expressions  for  3  e  and  5  e  follow  from  the  Prandtl-  Reuss  flow  equa- 
x  x  x  y 


» v  1  )-> 


■»»  mm 
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tions  (1-3).  By  considering  ^(e^-e  )  we  find 


B  v  -  ^  u  =  (2X/v  )  a  -  (l/>»)  *  c 
yx  xx  x  •  '  *  _ 

Similarly  the  relation  for  c  e  yields 

x  xy 


(3.9) 


3  v  +  5  u  =  -  (2\/v  )  a  +  (1/n)  ^  c 
xx  yx  =  xy  x  xy 

For  plane  strain  the  Tresca  yield  condition  can  be  expressed  in  the 

form 

2  2  ,2 
a  +  a  =  k 
xy 

Finally,  the  boundary  conditions  at  y  =  0  are 


x  <  0 
x  2  0 


cy  =  °  , 
v  =  0  , 


a  =0 
xy 


xy 


=  0 


(3.10) 


(3.11) 

(3.12) 

(3.13) 


It  is  noted  that  the  yield  condition  (3.11)  is  satisfied  identically 


by 


c  -  -  k  cosiu 


<J  =  -  k  siniu 
xy 


(3.14) 

(3.15) 

By  introducing  these  expressions  in  (3.7)  -  (3.10),  and  subsequently 
eliminating  (X/v  )  we  arrive  at  a  system  of  equations  of  the  general  form 


4  3u .  3u . 

I  L.  .  r-1  +  -rr  -  °  , 

j-i.  *  & 


(3.16) 


where 


ux  =  w, 


U2  = 


U3  =  Ux 


U.  =  V  , 

4  x 


(3.17) 


and  the  matrix  L  with  components  L_  is  defined  by 


I 

_ I 
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2k(  l+v)cosx 


k  cosx 


2k  ( 1+v) 
3  cos* 


-  tanx 


l+v)M^tanuj  -■§•  u(l+v)M^ 


k  _ 1_* 

U  cosx 


3 ( l-2v ) 
2u(l+v) 


-  2  tan«y 


(3.18) 


0  -  1  0 
L  0  2u(l+v)  1  U  J 

Equation  (3.16)  is  of  the  general  form  stated  by  Eq.(A.l)  of  Appendix 
A.  Just  as  for  the  anti-plane  case  we  seek  simple  wave  solutions  of  Eq. 
(3.16).  Unfortunately,  for  the  plane  strain  case  the  system  is  too  compli¬ 
cated  to  yield  an  analytical  solution.  It  is,  however,  possible  to  obtain 
simple  wave  solutions  for  small  values  of  M,  where  M  is  defined  by  (1.10), 
by  using  a  perturbation  method  which  is  discussed  in  some  detail  in 
Appendix  B. 

Expressions  for  the  characteristic  directions  v  *  dx/dy  at  small 
values  of  M  are  given  by  (B.17)  and  (B.18).  The  relevant  characteristic 
direction  for  y  2  0  is 


v(l)  _  1-sinx  _  M  1-v 
3  cosw  vl+sinu) 


V  2 
■)  +  0(0 


(3.19) 


where  tu  depends  analytically  on  M. 


As  discussed  in  Appendix  B,  the  following  relations  on  the  characteristic 
curves  defined  by  are  obtained: 


2  2 
C  =  -  yk  (1+v)  <u  +  +  0(M  ) 


(3.20) 


ux  -  -  ^  -q(l-v)  (l+sinuj  )^  -  +  u°  +  O(M^) 


(3.21) 
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v  =  %  (l-v)^(l-sinsu  +v°  +  0(M2) 
x  nM  x 


(3.22) 


Substitution  of  (3 .20)- (3 . 22)  into  the  first  of  (3.16)  yields 


1-sirsu 


1-v  yi 

■sinaj  /  J 


5  m  +  3  •!)  =0 

x  y 


(3.23) 


Along  the  family  of  characteristic  curves  defined  by  ,  (3.23) 
reduces  to  duu  =  0  ,  i.e.,  w  =  constant..  Consequently  the  relation  (3.19) 
given  by  =  dx/dy  ,  can  be  integrated  to  yield 


• M  (i&r)  ] * *  °<e> 


(3.24) 


where  3>(S)  is  an  arbitrary  function  of  the  polar  angle  9.  For  the  appli¬ 
cation  at  hand  the  appropriate  choice  is  t>(6)  s  0  for  a  central  field  of 
characteristics . 

Equations  (3.24)  and  12.21)  now  yield  the  following  relations 


-k  2 

sinsu  =  -  cos20  -  M  i_2(l-v)j*  sine  sin26  +  0(M  ) 

r  2 
cosu)  =  sin26  -  M  L2(l-v)]  sin0  cos26  +  0(M  ) 


(3.25) 


(3.26) 


and  hence 


u  =  20  -  -j  -  M  [2(l-v)]^  sin6  +  0(M2) 


(3.27) 


Substitution  of  this  result  into  (3.14)  and  (3.15)  yields 

a  =  k  fcos20  +  M  [2(l-v)l^  sin6  sin26} 
xy 


(3.28) 


o  =  k  {-sin29  +  M  L2(l-v)X*  sin0  cos2©} 


(3.29) 


The  average  pressure  follows  from  (3.21)  as 


c  =  Oq  -  |k(l-K>)  {29  -  ^  .  M  [2(l-v)p  sin0 }  +  0(M2) 


(3.30) 


i 
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By  combining  (3.2),  (3.3)  and  (3.4)  we  find 


2(1+0 


3c  \  j,  o 

,  t  +  k  -r  -26-sin28  ,  +  Mk  _2(l-v)_  cosS  sin28  +  0(M  ) 
2(L+v)  1 


(3.31) 


y  2  ( l+v ) 


.  -  °— -  +  k  (?  -26+sin26  + 2Mk  [2(l-v)]^  sin^S  +  0(M2) 

L  (  1+V )  >  2.  / 


(3.32) 


Cq  -2vk  [26  -  |  -  M[2(l-v)]^  sinG]  +  0(M2) 


(3.33) 


In  these  expressions  the  constant  is  still  undetermined.  Additional 
relations  can,  however,  be  obtained  from  the  boundary  conditions  (3.12)  and 
(3.13).  Clearly  (3.28)  and  (3.32)  can  neither  satisfy  the  conditions  stated 
by  (3.12)  nor  by  (3.13)  which  suggests  that  (3.28)  -  (3.33)  are  valid  only 

if  "k  k 

in  a  domain  3^  s  3  4  9^  ,  while  other  solutions  hold  for  0  £  8  £  3^  and 

k 

$2  ^  8  £  n  .  Appropriate  solution  foT  the  latter  domains  are  again  constant 


*  * 

For  both  0  ^  9  £  and  0O  £  9  £  rr  constant  states  imply  c  *  0  . 

12  r  xy 

k 

Equation  (3.28)  and  the  condition  that  J  should  be  continuous  at  6  =  8, 
^  xy  1 

★ 

and  8  =  then  yields  the  equation 

cos28  +  M  L2(l-v) sin8  sin26  =  0 


(3.34) 


which  to  first  order  in  M  has  the  solutions 


9*  =  J  +  -|(l-v>^  M  +  0(M2) 


(3.35) 


k  3tt  1  %  2 

92  =  T  +  2  (1"V)  M  +  °(M  } 


(3.36) 
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Next,  consider  (3.16)  for  i  =  3,  and  eliminate  tc/rx.tu  rx  and  fv  >x  bv 

x  x 

using  (3.20)  -  (3.22).  The  result  is 

=  f(«)  $  (3.47) 

by  vox.  tx 

where  f(»)  can  easily  be  determined.  Equation  (3.47)  can  now  be  integrated 
with  respect  to  x,  and  w  can  subsequently  be  eliminated  by  using  (3.27). 

The  result  of  these  manipulations  is 

k  2  k  2  1  2k  r-  _  ft 

u^  =  -  —  £n(2sin  8)  -  2— (1-v)  cos  8  +  —  —  L2(l-v)i^  Lcoso  +  £n(tan  j 


+  2ur  (y) 

Here  HO)  =  0,  but  otherwise  Hy)  is  an  arbitrary  function.  The  strain 
component  «  then  follows  from  (3.44)  and  (3.47)  as 


(3.48) 


xy 


-  -(1-v)  cos28  -  2n(2sin28)  +  £  -  [2(l-v)]*  [2cos6  +  In  (tan  |)] 


+  i  ux  +  y(y) 

At  8=8^,  Eq.(3.44)  yields 

vx  =  vx  +  M  +  0(M) 

*  ★ 
Since  v  Is  continuous  at  8  =  8  ,  and  uniform  in  the  domain  0  £  6  £  8,  , 
x  1 

and  since  by  symmetry  considerations  vx  must  vanish  at  6  =  0,  it  follows 
that 


(3.49) 


(3.50) 


vx  “  ’  M  ^(1_V)%  +  °<M> 


(3.51) 


By  evaluating  (3.43),  (3.46)  and  (3.49)  at  8  =  8^  we  then  obtain  for 
0  *  0  *  &*: 


o  1  2k,. 
t  *  u  -  —  — (1-v) 
x  x  M  n 


(3.52) 


f 
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t  -  *  \  +  ~(1-2'j) (1+r)  +  I  — (l-v)* 
y  x  h,  M  > 


(3.53) 


xy 


=  7(y) 


(3.54' 


A  computation  of  u^  can  be  based  on  the  assumption  that  at  s  =  0 

the  domain  of  validity  of  (3.52)  -  (3.54)  borders  directly  on  the  zone  of 

elastic  deformation,  and  that  a  ,  a  and  a  are  continuous  at  the  elastic- 

x  y  x 

plastic  boundary  .  If  in  the  plastic  zone  is  represented  by  (3.40),  and 
outside  the  plastic  zone  is  computed  by  using  Hooke's  law  with  s^  and  a 
as  defined  by  (3.52)  and  (3.53),  then  continuity  of  yields 


<  -  3  T<l-V)%  +  z  [i  - 


(3.55) 


The  strains  in  the  domain  0  £  3  S  then  reduced  to 


e  =  -  Q  -  V(l-HT)] 
X  p.  z 


(3.56) 


«  -  ^  C-f  +  (l-^)d-KT)] 


(3.57) 


By  evaluating  (3.43),  (3.46)  and  (3.49)  at  0  =  6^  ,  and  using  (3.51) 
and  (3.55)  we  find  for  9^  £  0  £  n  : 

«x  -  *  Cf  -  1  +  d-v//2] 
ay  =  ^  i-j  +  (l-v)(4-/2)/2] 


(3.58) 

(3.59) 


1  Jr  ]>•  %  V 

*  -  -  =  -  *d-v)(l-l//2)  +  -  L2(  1-v)  J*  Zn(tan3n/8)  +  f  (y)  (3.60) 

xy  “  U  li 


The  crack  opening  angle,  which  is  defined  by  tana  Hvx|»  follows  from  (3.44)  by 

it 

substituting  0  =  The  result  is 


1  4k., 

tana  =  —  — ( 1-v) 

M  p, 


(3.61) 
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The  expressions  for  the  stresses  given  by  (3 . 31)- (3 . 33)  show  only  a 
weak  dependence  on  M.  In  the  limit  of  vanishing  crack- tip  speed  (M— 0) 
these  expressions  reduce  to  the  corresponding  quasi-static  results  as, 
for  example,  presented  by  Rice  and  Tracey  (1973).  In  the  range 

k  k 

<  S  <  $2  *  t*le  strains,  given  by  (3.43),  (3.46)  and  (3.49)  strongly 
depend  on  M  due  to  the  presence  of  terms  of  order  0(M  ^) .  The  strain 
and  e ^  have  been  plotted  in  Figs.  5  and  6  for  M  *  0.1  and  M  =  0.01. 

In  the  limit  M  —  0  the  strains  become  unbounded. 

The  results  presented  in  this  section  do  not  show  elastic  unloading. 

For  the  quasi-static  problem  it  has  been  shown  by  Rice,  Drugan  and  Sham 
(1979)  that  an  elastic  unloading  zone  must  exist  in  the  immediate  vicinity 
of  the  crack  tip. 

Thus,  the  near-tip  in-plane  asymptotic  results  presented  in  this 
Section,  which  take  into  account  dynamic  effects,  display  the  same  anomalies 
that  were  observed  in  Section  2  for  the  case  of  anti-plane  strain.  There 
is  no  uniform  transition  from  the  dynamic  to  the  quasi- static  results  as  the 
crack- tip  speed  becomes  very  small,  and  the  dynamic  solutions  do  not  show 
elastic  unloading  in  the  immediate  vicinity  of  the  crack  tip.  Just  as  for 
the  Mode-Ill  case  it  is  speculated  that  the  actual  transition  from  dynamic 
to  quasi-static  results  is  achieved  because  the  zone  of  validity  of  the 
results  of  this  section  shrinks  on  the  crack  tip  in  the  limit  M  -  0. 
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APPENDIX  A 


Simple  Nave  Solutions 

The  systems  of  equations  chat  appear  in  this  paper  are  of  the  general 


-U  .  pU^ 

Lij  rx  By  ° 


i  =  1,2. 


(A.  1) 


Equation  (A.l)  defines  a  system  of  n  simultaneous  first-order  partial 

differential  equations  for  the  dependent  variables  u, ,u,,  -  -  -  -  ,  u  . 

1  n 

The  equations  are  homogeneous,  and  are  functions  of  x,y  and  u;  ,  but 
not  of  the  derivatives  of  the  u^.  Such  equations  are  called  quas i- linear , 

As  discussed  by  e.g.  Bland  (1969)  solutions  of  Eq.(A.l)  are  simple  wave 
solutions  if  all  the  u^  are  functions  of  one  another,  or  equivalently, 
if  all  u.  are  functions  of  another  variable,  sav  3  ,  For  u.  =  u.(3)  we 

l  J  li 


*u.  -  Au.  .  du. 

i  /  B*P  1  /  A©  ,  r 

— r—  ■  u.  —  ,  — r~  =  u.  y—  ,  where  u.  =  — — 

Bx  i  Bx  py  i  By  i  dec 


(A. 2) 


Equation  (A.l)  can  then  be  rewritten  as 


*  /  Bed  ,  ,  Bcd  n 

Z  L.  .u.  —  +  u.  “  =  0 

ij  j/  i  By 


(A. 3) 


Along  any  curve  in  the  xy- plane  we  have 


dx  ■“  +  dy  ~  k  dco 

Bx  J  By 


(A. 4) 


Hence  Eq.(A.3)  define  total  differentials  along  curves  defined  by 


n  n 

I  L.  u  '  L  L.  .du 

1-1  J  _  1-1  J 


dx  m  1=1 
dy  u 


Equations  (A. 3)  are  true  for  all  values  of  i,  and  therefore,  the  right- 


hand  side  of  Eq.(A.5)  is  the  same  for  all  values  of  i.  Hence 


j  =  l 


(L.  .  -  V  5.  .)  du.  =  0 

7  1  1  1  '  1 


J-J 


J 


i  -  1.  2,  . 


»  n 


(A.  6) 


This  system  of  equations  has  non-trivial  solutions  only  if  the  determinant 
of  the  coefficients  vanishes,  i.e.,  if 


L.  .  -  v  6.  .  =0 

ij  ij  1 

Equations  (A. 7)  is  recognized  as  the  equation  for  the  slopes  V  =  dx/dy 
of  the  characteristics. 

The  relations  holding  on  a  characteristic  of  a  simple  wave  follow 
from  Eqs.(A.6).  Since  the  du .  are  the  right  eigenvectors  of  L. .  - 


(A. 7) 


J 


ij 


'ij' 


they  are  proportional  to  the  cofactors  of  any  row  of  L.  -  V  6 

ij  ij 

provided  that  such  cofactors  are  not  all  zero.  Thus  we  have 

du.  du.  du 

2  n 


(A. 8) 


where  C.  ,  i  =  1,2 

i  ’ 


,  n  ,  are  such  cofactors. 


\ 
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APPEND  LX  B 
Perturbation  Method 


The  system  of  governing  equations  for  the  case  of  plane  strain,  given 
by  (3.16),  is  of  the  general  form  (A.l).  We  again  seek  simple  wave 
solutions.  Formally  the  equation  for  the  characteristic  curves  is  given 
by  Eq.(A.7).  This  equation  is,  however,  much  too  complicated  to  be  solved 
analytically.  Thus,  it  was  attempted  to  obtain  simple  wave  solutions  for 
small  values  of  M,  where  M  is  defined  by  (1.20). 

For  M  =  0  ,  the  fourth-order  equation  given  by  (A. 7),  with  defined 
by  (3.18)  yields  the  following  double  roots: 


o  1+s  inu) 

v ,  =  -  - g— 

1  COSO 


(B.l) 


o  =  l-sinu)° 
'2  cosoi0 


(B.2) 


f 
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Here  the  superscripts  zero  indicate  that  M  =  0  .  The  multiplicity  of 
the  roots  necessitates  some  special  provisions  before  a  perturbation 
procedure  can  be  attempted.  If  expansions  for  small  M  were  directly 
substituted  in  (A. 7)  the  result  would  be  a  perturbation  of  (B.l)  and 
(B.2),  while  the  two  other  different  roots  would  not  be  obtained.  This 
difficulty  can  be  circumvented  by  preceding  the  perturbation  procedure 
by  a  transformation  which  would  reduce  L„  at  M  *  0  to  the  following 
canonical  or  Jordan  form- 


1+sinw 


0 


l+sinw 


(B.3) 


1-sinu) 


1-sinaj 


For  the  case  of  multiple  eigenvalues  formal  methods  are  available  to 
reduce  a  matrix  to  its  canonical  form.  Here  it  suffices  to  state  the 
transition  to  the  Jordan  representation  can  be  achieved  by  introducing 
the  new  variables  i,  a,  ux  and  vx  by  the  following  substitutions 


t  dS  +  a  u_  di 

k  L-v  k  1-v  x 


(B.4) 


,  2  1+v  .*  2  1-fv  . 

da  *  -  —  u  - —  do  -  T  u  - —  cv 

3  l-v  3  *  l-v  x 


(B.5) 


1+sinai  , ~  .  v  .  l-simu  .  v 

-  duu  +  - —  da  + -  du  +  - —  dv 

cosut  l-v  cosa>  x  l-v  x 


(B.  6) 
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dv  =  duj  +  dv 

X  X 


(B .  7) 


For  the  purposes  of  the  perturbation  procedure  we  introduce  in  addition 


dv  »  i  dv  ,  da  =  7  dc 
M  x  ’  M 


(B.8) 


With  the  substitutions  (B.4)  -  (B.8)  the  system  of  equations  corresponding 
to  (A. 6)  becomes 


-  [v  +  (1  -  |  M2)  +  M(1  -  £  sin«)  da  + 


2  3 

.  v  ..2  1-sinu)  .  -  v  M 

+  x  M  -  du  -  - —  t—  sinuu  dv  =  0 

2  cosuj  x  1-v  2  x 


(B.9) 


M  di  -  fv  +  (1  +  |  M2)  i±aiS2il  dS  -  £  M2  £  =  0 

l-sirau  L  2  cosuj  j  2  cosuj  x 


COSUJ  X 


(B.10) 


2  3 

v  ..2  1+sinuu  ,  -  .  v  M  . 

-  7  M  -  duj  +  - - r  sinuj  da  + 

2  cosuj  1-v  2 


2  ..2 


+  f  "V  +  VosuT  (1  '  I  m2)]  dux  +  M(1  +  T  sin“)  dvx  =  0  (B- 11) 


v  m2  l^sinu,  d~  +  M  d~  +  f  -  y  +  IzliSS.  (i  +  v  m2 

2  cosuj  1+sinuj  x  L  cosuj  v  2 


d+7  M2)]  dv^  =  0  (B.  12) 


The  canonical  form  given  by  (B.3)  is  recognized  as  appearing  in  the 
limit  M  =  0. 

At  this  stage  we  can  consider  expansions  in  terms  of  powers  of  M: 


V°  +  M 


(B. 13) 


Lij  "  Lij  +  M 

together  with  analogous  expansions  for  duj,  da,  du^  and  dv^,  while  uj  is 
assumed  to  depend  analytically  on  M. 


( B .  14 ) 
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On  the  basis  of  Eq.(A.7)  we  obtain  in  the  zero'th  approximation  the 
results  given  by  (B.l)  and  (B.2).  Considerations  based  on  (A. 8)  yield 
the  following  relations  along  the  characteristics 


,  1-sirwj 

du  =  ±  -  dv  , 

x  ox’ 

COSU) 


d»  »  dc  *  0 


(B . 15) 


The  first  order  approximations  are  obtained  as 


y(D  _  _  l+sinq;  ±  M  \  n(M2 

1,2  ‘  >l-sin“«/  ' 


(B.16) 


,(1)  _  1-sing;  /  1-v  y  2 

3.4  - ™  17*15;  +°(M) 


(B. 17) 


The  eigenvectors  corresponding  to  Y^  ’  lead  to  the  following  relations 


=  dv  =0  ,  dc  =  (-r^— r — )  di 

x  ’  \l-sinu)/ 


(B.18  ) 


In  terms  of  the  original  variables  we  find 


.  3  dc 

=  2k  ISJ 


(B.19  ) 


i  =  [  (1-v)*  (1+sinw)^  -  —1  dm 

X  l  n 


(B.2Q) 


dvx  -  -  ^  (1-v)*  (1-sirw)*  dm 

These  expressions  can  easily  be  integrated.  In  the  following  we  list 
not  only  the  expressions  corresponding  to  Y^  (upper  signs),  but  also 
the  ones  corresponding  to  Y^^  (lower  signs); 


(B.21) 


u  *  ?  ~  (1-v)^  (1  ^  sinui)^  ?  —  +  u° 

X  liM  j.  X 


(B.22) 


* 


(B . 23  ) 


c  =  =  7  k  (1+v)  «  +  a 

J  A 


vx  =  "  (1'V)^  (1  x  sin*)J5  +  \ 

where  u°  ,  a  and  v°  are  arbitrary  constants, 
x  o  x 

It  is  of  interest  to  determine  \  from  the  relation 


(B.24) 


2\kr - t —  (-  sin8  +  cos8  3  u  'j 

sinui  \  9  x  9  x  / 

Corresponding  to  Y^  and  Y^1"*  we  find 
Xr  -  -  <  0 

kiM 

Analogously,  we  find  corresponding  to  Y^(Upper  sign)  and 
(lower  sign) 


(B.25 ) 


(B.26) 


CT  =  ±yk(l+v)*+ao 


ux  =  ~  ~  ^  (1  *  v)^  (1?  sinm;^  T  v  »j  +  u 


(B.  27) 
(3.28) 


vv  =  ^  <l-v)*  (1  ±  sinui)^  +  v° 
X  x 


(B.29) 


[2(l-v)]^  . 

\r - — <-*—  >  o 

kiM 


(B.30) 


